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Abstract—A proper set of models are carefully selected from 
the literatures to form a wholistic perspective in simulating the 
wave propagation in evaporation ducts, under a close-to-realistic 
environment. The split-step Fourier (SSF) propagation algorithm 
is applied to analyze the propagation properties in evaporation 
ducts above a rough sea surface, with turbulence in the ducts. 
The scattering effect of the rough sea surface is modeled with 
an effective reflection coefficient, including the Miller-Brown 
correction factor for surface roughness. The turbulence effect is 
simulated with a two-dimensional Kolmogorov power spectrum 
of refractive index fluctuation, which is derived from the three- 
dimensional Kolmogorov power spectrum using the Wiener- 
Khinchin theorem, including the anisotropic effect of the outer 
scales. The average /-profile, the outer scales of turbulence, and 
the structure constant of the refractive index are categorized 
under different atmospheric conditions. Simulation results are 
presented and discussed over all possible atmospheric conditions 
to better understand the effect of turbulence on the evaporation 
ducts. 


Index Terms—Propagation, split-step Fourier (SSF) algorithm, 
evaporation duct, rough surface, turbulence, Kolmogorov power 
spectrum, refractive index fluctuation. 


I. INTRODUCTION 


Measurement of centimeter-wave propagation indicates that 
the difference between the field values inside and above an 
evaporation duct is much smaller than the prediction using 
a simple model [1]. Such difference may be explained as 
follows [1], [2]: (1) The height and strength of an evaporation 
duct can vary with distance, and the wave within an irregular 
duct is less focused. (2) Roughness of sea surface redirects a 
large portion of wave energy into the space above the duct. 
(3) Turbulent fluctuations in the atmosphere tends to broaden 
the angular spectrum of waves and increase leakage from the 
duct. Scattering by turbulent atmospheric fluctuations is the 
strongest mechanism, at least in certain cases. 

For surface-based radars, non-standard tropospheric refrac- 
tion such as ducting may result in an anomalous propagation 
condition, bending the radar ray from the anticipated direction 
[3]. Anomalous propagation echoes have been observed in fair 
weather conditions; more frequently in summer than in winter, 
more frequently in tropical areas than in the middle and higher 
latitudes [4]. 

Ducting may seriously affect radio communications links, 
extend the radar range of a target at or near the sea surface 
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to be considerably longer than if in the free-space [5]. Evap- 
oration duct is a surface duct appearing over water bodies, 
featuring rapid decrease of humidity with height. It is the most 
common type of anomalous propagation over oceans and other 
large bodies of water. The worldwide mean of the evaporation 
duct height (EDH) is about 10 m [6], sometimes reaches 30 m 
or 40 m [7]. The EDH increases during the summer months 
and during the daytime [8]. With the evaporation ducts so 
close to the sea surface, the latter has very strong effects on 
the signal propagating in the ducts. 

The effects of duct on communications links have been 
studied [9], [10]. Possible effects of tropospheric ducting 
on short-distance wireless communications have also been 
reported [11]. The effects of evaporation and surface-based 
ducts on microwave path loss have been studied [12]. In [13], 
the path-loss variations are compared with the monthly and 
seasonal mean parameters of surface ducts, based on a two- 
year statistical study of surface-duct formation at Istanbul, over 
a WCDMA-FDD downlink band [14]. Ducts frequently appear 
in the coastal areas where the horizontal change of refractivity 
can not be neglected. In [12], a comparison is made between 
a range independent surface-based duct and a mixed land-sea 
path, using two M -profiles. 

Propagation experiments have been conducted in the 
Mediterranean and the North Sea areas [15], [16]. The received 
signal is stronger than expected, which is a sign that the 
evaporation ducting effect may dominate over the gaseous 
absorption in the 10-20 GHz band. Measurements at 18 GHz 
have been conducted to demonstrate the enhancement on radar 
coverage and range detection due to an evaporation duct [17], 
[7]. While the predictions and the measurements seem to agree 
reasonably well, the predictions tend to underestimate the path 
loss when the duct heights are low. Although the prediction 
model incorporates the Kirchhoff theory to account for the 
scattering from a rough sea surface, it does not consider the 
scattering due to the turbulent fluctuations in the refractive 
index [18]. 

An evaporation duct can significantly increase the signal 
strength near the sea surface for over-the-horizon paths, at 
frequencies above 3 GHz. As the frequency is higher than 
10 GHz, the signal strength is reduced by scattering from 
the wind-driven rough sea surface. Such signal reduction 
is more obvious under strong ducting conditions and large 
sea-wave heights [19]. In [19], an MLAYER propagation 
model, including the Miller-Brown reflection reduction factor 
(RRF), is used to compare with the measurements. In all the 
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cases studied in [19], the model including rough sea surface 
provides a closer match with the measurement results than 
the model with smooth sea surface, but the former model still 
underestimates the path loss by 3 to 12 dB. The refractive- 
index fluctuation may partly account for this underestimate. 

The effects of air turbulence have been analyzed by applying 
a small-perturbation method to a mode-based model [18], or 
applying a phase-screen method to a split-step Fourier (SSF) 
model [20], [21]. Different models have been compared on 
a series of over-water measurements, in a nearly standard at- 
mosphere [21]. The model including the turbulence scattering 
performs better than the other models, in the sense of matching 
with the measured data. At higher frequencies, the effects of 
turbulent scattering become even stronger. 

In [22], two approaches are proposed to model the surface- 
layer turbulence effects in a marine environment at microwave 
frequencies. One is the SSF method with phase screen, the 
other is a direct computation of the log-amplitude variance of 
signals. Their results are compared with the observation, and 
the effect of the power-law constant on the first approach is 
also discussed. 

In [23], the phase screen method is adopted in the split- 
step Fourier algorithm to simulate the turbulence effect in an 
evaporation duct. It is assumed that the atmosphere is under the 
neutral condition, the turbulence follows the one-dimensional 
Kolmogorov spectrum, the structure constant and the outer 
scale are constant, and the sea surface is smooth [23]. The 
simulation results show that the refractive-index fluctuation 
leads to energy leakage away from the duct, causing an 
additional path loss within the duct. The height distribution of 
radiated power becomes more uniform than that without the 
fluctuations. The power leakage along the evaporation duct is 
not uniform. It is also observed that the received power level 
above the duct fluctuates more significantly than within the 
duct when the refractive-index fluctuation is considered. 

Diffraction by a relatively smooth sea surface can be de- 
scribed by introducing an effective reflection coefficient [24] 
into the boundary condition, which is the flat-surface Fresnel 
reflection coefficient multiplied by a roughness reduction 
factor (RRF). In general, a rough sea surface tends to destroy 
the trapping property of the duct structure and changes the 
path-loss pattern. In [25], the effects of wind direction on 
electromagnetic waves propagation over rough sea surface are 
studied by mixing two approaches applicable at large and 
small roughness, respectively. 

The interaction among sea waves is nonlinear, rendering 
the sea-surface statistics different from a Gaussian distribution 
[26]. It has been demonstrated that the shadowing effect 
prevails over the non-Gaussian statistics of a sea surface, 
especially in the presence of an evaporation duct at the 
centimeter wavelength [27]. 

The wave equation can be reduced to a parabolic equation 
(PE), which is then solved using numerical techniques like the 
split-step Fourier (SSF) algorithm [28], finite difference (FD) 
algorithm [24], and finite element (FE) algorithm [29]. The 
FD and FE algorithms are more flexible to implement various 
boundary conditions. The SSF algorithm is numerically stable 
and allows a larger step size to compute the field distribution 


at long distances. 

Reasonable agreement between measured and predicted 
field strength using the PE model has been achieved in a 
tropical maritime environment [30]. It is shown that the PE 
model gives a good estimation of path factor over a wide 
range of frequencies (X, Ka and W bands) and a variety of 
atmospheric conditions [31]. The difficulties in predicting the 
field strength of maritime links and the importance of higher 
M -inversion can be found in [32]. 

Numerical analysis indicates that the microwave inside an 
evaporation duct suffers an additional path loss in the presence 
of turbulence [33]. Monte-Carlo method has been used to 
study the scattering properties of scalar waves in randomly 
fluctuating slabs with an exponential spatial correlation, and it 
has also been extended to random media with non-exponential 
spatial correlation and background inhomogeneity [34], [35]. 

The prediction accuracy of the PE models is constrained 
by the environmental parameters, especially in littoral envi- 
ronments [12], [36]. Three-dimensional PE models have been 
developed [37], [38], covering the sea-surface roughness effect 
[25] and the evaporation ducts [39]. PE-based time-dependent 
propagation model has also been explored [40]. 

Although many of the key factors on evaporation ducts, 
like refractive-index profile, rough surface, wind profile, tur- 
bulence, have been analyzed or modeled; few literatures, if not 
none, integrate all the relevant factors together to analyze or 
simulate wave propagation in evaporation ducts under a close- 
to-realistic environment. In this work, we carefully review 
many existing models about these factors, and choose a proper 
set of models to simulate the wave propagation in evaporation 
ducts, under all possible atmospheric conditions. During the 
review process, it is observed that the wind profile, which can 
be categorized with a stability index of the atmosphere, seems 
to be a common factor behind all these models. As a result, this 
work provides a wholistic perspective on the subject issues. 
The simulation results may also help reveal more values of 
the existing models and inspire more ideas from the readers. 

In this work, the SSF algorithm is used to simulate wave 
propagation over long distances. The Miuller-Brown [42] for- 
mula of RRF is adopted to model the rough sea surface. 
Monte-Carlo method is used to generate profiles of refractive- 
index fluctuations to simulate the turbulence in the atmosphere. 
The relevant models are described in Section I, the proper 
ranges of parameters in these models are estimated in Section 
III, the simulation of practical evaporation ducts is presented 
and discussed in Section IV, followed by the conclusion. 


II. CONSTRUCTION OF MODELS 
A. Propagation Model 


The parabolic equation is a paraxial approximation to the 
Helmholtz wave equation, under the assumption that the elec- 
tromagnetic wave predominantly propagates in the horizon- 
tal direction, 2. The two-dimensional narrow-angle forward- 
scattering scalar form of the parabolic equation takes the form 
[28] 


8? 
1] + aa} u(x,z) (1) 


z = [m?(z, z) — 
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where k is the wavenumber in the background medium, m = 
1+.M x10~°, with the modified refractivity, M, related to the 
refractivity, N, as M = N+0.157z. Eqn. (1) is very accurate 
with the propagation angle about the horizontal direction falls 
within |0| < 15° [24]. 

The split-step solution to (1) can be expressed as [28] 


u(x + Az, z) = e71 (k/2 (m -1)Az 
pe Faeee))) (2) 


where p = ksin@ is the vertical wavenumber. The solution 
u(x, z) and its spectrum U(x,p) are related by the Fourier 
transform F{-} and its inverse F~'{-} as 


iS -= l. u(x, zjedz 


“val” 


To represent wave propagation over a perfect electric con- 
ductor (PEC) surface, either a Dirichlet (horizontal polariza- 
tion) or a Neumann (vertical polarization) boundary condition 
is imposed. To represent wave propagation over a smooth 
impedance surface, the boundary condition can be specified 
as [43] 


Fultz) 
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where x = 90° — 9; is the grazing angle, Reg is the Fresnel’s 


reflection coefficient, R} or RY, with 


sin X — y Esr — COS* X cos? y 
Ia r. (5) 
Sin X + 4/ €sy — COS 
2 
Esr SIN X — \/ Esr — COS 
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for the horizontal and the vertical polarization, respectively; 
Esr = €s/€9 = €r — J600A 


is the effective dielectric coefficient of the sea sea surface, o 
and e, are the conductivity and relative permittivity, respec- 
tively, of the sea water [44]. 

The boundary condition can be incorporated into the stan- 
dard split-step PE approach through the use of mixed Fourier 
transform [28], [43]. The path loss (PL) is defined as [24] 


dya 7 dB (7) 


PL = 20 log sary 


where A is the wavelength of the wave, and x is the propaga- 
tion range. 


B. Effective Reflection Coefficient of Rough Sea Surface 


The roughness effect can be incorporated by multiplying the 
Fresnel’s reflection coefficient, Rr, by a roughness reduction 
factor (RRF), pmr, to obtain an effective reflection coefficient 
[24] 


Ref = PMBEF 


where 
pup = e7? /? Ig (97/2) (8) 


is the Miller-Brown roughness reduction factor [42], Io(a@) is 
the zeroth-order modified Bessel function of the first kind. The 
relevant parameters have the form 


y = 2kop sin x 
onr=bIi x0 Ua (9) 


where øp is the standard deviation of the sea wave, k is the 
wavenumber of the signal, and Uj is the reference wind speed 
at 10 m above the sea level. Note that the RRF in (8) affects 
the magnitude of the complex Fresnel reflection coefficient, 
but does not account for the shadowing effect [45] or the 
geometry of the sea surface [25]. At near grazing incidence, 
the shadowing effect may have stronger effect on the wave 
propagation [45]. 


C. Refractive Index Fluctuation 


The total refractive-index is expressed as 
nz) + nsz, z) (10) 


where 7(z) is the average refractive index, which is a function 
of altitude, and n(x, z) is the fluctuation of the refractive- 
index. The turbulence is approximated as frozen [46] during 
the period the electromagnetic wave propagates from the 
transmitting antenna to the receiving one. By applying the 
Monte-Carlo technique, one sample profile of refractive-index 
fluctuation, of a two-dimensional turbulence, can be realized 
as [34], [41] 


Plea) = DoS 


p=1q=1 
sin(Cept + Czqgz) + pis) COS (Cent + zq2) (11) 


where s indicates the sth realization, ¢;, and Cz, are the pth 
and qth discretized wavenumber, respectively, a&pq and bpg are 


random numbers with variance a which can be expressed 
as 


i t.2) = 


= 4ACrpAGzq Fin (Gap, Szq) (12) 


where F7, (Cz, bite is the two-dimensional power spectral den- 
sity of the refractive-index fluctuation. 

The three-dimensional power spectral density of refractive- 
index in an isotropic Kolmogorov turbulence is [47] 


0.03302 
(Grate tli 


Where Lo is the outer scale, and C2 is the structure constant. 
By the Wiener-Khinchin theorem [48], the two-dimensional 


oc) (13) 
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power spectral density, independent of the y coordinate, can 
be derived as 


FE (n, G) = J BE (Ca, Cys Cde, 


0.055507 
= ae n (14) 
(Gres li 
which is an isotropic Kolmogorov turbulence. Next, follow 
the same approach as in [47] to modify (14) to be a two- 
dimensional anisotropic Kolmogorov spectrum 


0:0555C7(Logling 
(C2Lon + C2L6, + 1)4/8 


where Lo, and Lop are the outer scales in the vertical and the 
horizontal directions, respectively. 


FE (e, G) = (15) 


III. ESTIMATION OF PARAMETERS 
A. Average Modified Refractivity 


The PJ model is adopted to generate profiles of modified 
refractivity under different atmospheric conditions as [49], [50] 


4(1 + z/2z0) 
(1 ar 162/Eno) 
a eo ee, 
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where Lmo is the Monin-Obukhov length, zą (n m) is the duct 
height, and zọ (in m) is the roughness length of the surface. 
Typically, zọ increases with the wind speed, ranging from 
1074 m on a calm open sea to 1078 m in coastal areas with 
off-sea wind [50]. The value of z/Lmo is used to categorize 
the atmospheric condition: Stable if z/Lmo > 0, unstable if 
2) Lines < 0 and neutral if z/Lire =Ù.: 


B. Monin-Obukhov Length 


The Monin-Obukhov length is essentially dependent upon 
the heat flux and the friction velocity, which can be treated 
as independent of height within the surface layer. The gen- 
eral characteristics of the atmosphere can be categorized in 
terms of z/Lmo as follows [50]: When Lmo is negative with 
very small magnitude ( z/Lmo becomes negative with large 
magnitude ), the heat convection dominates. When Lmo is 
negative with large magnitude ( z/Lmo becomes negative with 
small magnitude), the mechanical turbulence dominates. When 
z/Lmo = 0, only mechanical turbulence contributes. When 
z/Lmo is a small positive number, the mechanical turbulence 
is slightly damped by the temperature stratification. When 
z/Lmo is a large positive number, the mechanical turbulence 
is severely reduced by the temperature stratification. The 
number |z/Lmo| describes the relative significance between 
heat convection and mechanical turbulence during the daytime, 
and to what extent the stratification suppresses the mechanical 
turbulence during the night time, respectively. 


TABLE I 
PARAMETERS DERIVED FROM MEASUREMENTS [52]. 


Stability | Lmo (m) | ux (m/s) | zo m) zi (m) Exp. 
class no. 
62xio-> [- | 358 
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zi is the height of the boundary layer. 





The Monin-Obukhov length can be estimated from the 
Richardson number, R;, as [50] 
2) Linge = (17) 


The Richardson number fF; in the surface layer can be 
_ IJ Wd yI 


expressed as [50] 
9 w-7_ (,, 0.07 
— T (dup, / Oz)? B 


where g is the gravity acceleration, y = —OT'/0z is the lapse 
rate of temperature, ya = 9.8° C/km is the dry adiabatic lapse 
rate, vp 1s the mean wind velocity in the horizontal direction, 
and 


Ri (18) 


BE T> — Ti + 0.01 (22 = 21) 
2500(q2 — q1) 


is the estimated Bowen ratio [50], with q the specific humidity. 

In [51], seven stability conditions are defined and related 
to the Monin-Obukhov length: 10 < Lmo < 50 m for very 
stable (vs) condition, 50 < Lmo < 200 m for stable (s) 
condition, 200 < Lmo < 500 m for near-neutral stable 
(ns) condition, |Lmo|] > 500 m for neutral (n) condition, 
—500 < Lmo < —200 m for near-neutral unstable (nu) 
condition, —200 < Lmo < —100 m for unstable (u) condition, 
and —100 < Lmo < —50 m for very unstable (vu) condition. 

In [52], measurements were conducted in 2006, at site 
located 18 km from the west coast of Denmark in the North 
Sea. A set of mean parameters derived from the measurements 
are summarized in Table I, categorized under different atmo- 
spheric stability conditions just mentioned. 


(19) 


C. Outer Scale of Turbulence 


For turbulence in the free atmosphere, the outer scale Lo 
typically ranges from 10 to 100 m. The isotropic assumption 
is usually valid for well-developed turbulence far from any 
boundaries. In a surface layer, which lies in the lowest 10 
% of a planetary boundary layer [50], the horizontal outer 
scale might remain unchanged, but the vertical outer scale is 
usually reduced, rendering the turbulence no longer locally 
homogeneous and isotropic in this layer. The turbulence in 
the surface layer is dominated by the wind-shear. A first-order 
vertical outer scale in the surface layer can be expressed as 
[47], [53] 


Lin = katd,, "0.7" (20) 
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where k, is the von Karman’s constant, @,, 1s the normalized 
wind shear, and ¢, 1s the normalized turbulent kinetic energy 
dissipation rate. The von Karman’s constant is conventionally 
taken to be 0.4. Its value measured in the wind tunnel and 
the atmosphere ranges from 0.35 to 0.43 [50], [54]. The 
normalized wind shear is defined as [50] 


kaz Ou 
Om — aD. 


üp OZ 
and the normalized turbulent kinetic energy dissipation rate is 
defined as [50] 





(21) 





a= 2 (22) 


us 
where u is the wind velocity, u, is the friction velocity, and e€ is 
the dissipation rate. The profiles of m and ¢, under different 
conditions are [50] 


(Hoes) =. Oe a 
om = (1—152/Lmo)1/*, 2/Lmo < —2 

1 + 5z/Lmo;, Zf Lms > 0 

1 — 2) lanes 2) Line <0 
b= m (23) 

1 42.5 (z/Tmo)” |  2/Lino > 0 


The vertical outer scale of the turbulence in a surface layer, 
under different atmospheric conditions, are derived by substi- 
tuting (23) into (20): 


kaa =al] 


; —2 < z/Imo < 0 
(1 — 162z/Limo)?/® / 
B -ija 
Los = a, TET a (24) 
(1 —152/Lmo) 
yall Bafa) ip 


iraia T 
where ka = 0.4 [47]. 


D. Anisotropy of Turbulence 


Turbulence at large scales is usually anisotropic, with 
considerably different outer scales in the horizontal and the 
vertical directions. The anisotropy and the length of the outer 
scale is dependent upon the atmospheric stability [55], [56]. 
The anisotropy of the outer scale becomes more obvious 
with stronger temperature inversion and weaker wind, which 
correspond to a higher Richardson number and a more stably 
stratified atmosphere. In such a weather condition with stable 
air, the mechanical turbulence is strongly damped by strati- 
fication, hence the vertical outer scale become much smaller 
than the horizontal one. In [55], it is also confirmed that if the 
turbulent energy is mainly generated from mechanical origin, 
the eddies are small, the outer scale length tends to decrease 
with increased wind speed, and vice versa, which is much 
more obvious in the horizontal scales. 


The anisotropy has no significant influence on power fluc- 
tuation, but the angle-of-arrival [55]. In [50], the ratio of the 
horizontal to the vertical integral scale (similar to the outer 
scale) of wind fluctuations has been estimated. It is shown 
that with an increasing positive Richardson number, a much 
larger outer scale is expected in the horizontal direction. Under 
neutral condition (R; œ~ 0), the horizontal and the vertical 
outer scales are similar. In an unstable air (R; < 0), eddies 
tend to be vertically elongated. In [56], it is empirically found 
that Lon/Lov œ (30 + 10)R; under stable condition. The 
wide range of variation might be attributed to the fact that 
the anisotropy does not solely depend on Ri. 

The anisotropy of the turbulence can be related to the 
stability of atmosphere as [57] 


— 100 —2/ Limo) 
tz = D5 z/Lmo 
where a; = 4r Lon /Lov is the decay constant of the vertical 
coherence of wind components [50], [58]. Eqn. (25) implies 
10(1— z/Lmo) 
4r (0.5 — z/Lmo)’ 


z/Lmo < 0.5 (25) 


Lon/ Lov = 2) Lie, SO. (26) 


E. Structure Constant 


Turbulence is likely to become strong ( structure constant, 
C2, becomes large ) near the earth surface and in the clouds 
[59]. The magnitude of C? generally increases with increasing 
wavelength, because the contribution of moisture to C? signif- 
icantly increases with increasing wavelength. The magnitude 
of C? in near-millimeter-wave band can be larger than that in 
the infrared by orders of magnitude. The magnitude of C? in 
the visible and the infrared bands near the ground falls in the 
range from 10716 to 1071? m~?/3 [60]. 

The structure constant of the refractive-index fluctuation in 
the planetary boundary layer can be expressed as [47], [53] 


2 
2_ 274/3 O(n) 
Ca — a Lo, (Se ) 


where a? lies between 1.5 and 3.5 [47], [50]. Substituting (16) 
into (27), we have 


C2 = 2.814? x 1071? 


0.032 i ZdV 1 — 1624) Lane 


(27) 


(agg at] 
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z/Lmo < 0 (28) 


8(1 + 5z¢/Lmo) 
Cling 0 


where a? = 2.8 and k, = 0.4 are chosen [47]. 

In the experiments considered in [61], the moisture fluctu- 
ation plays a dominant role, as compared to temperature, in 
determining the microwave C2; and significantly affects the 
acoustic and optical C2 values, since the near-surface turbu- 
lence leads to relatively large moisture fluctuations. The over- 


land experiment shows a larger diurnal variation, while the 


0.082 , 24[5/Lmo + 1/(z + a E 
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moisture fluctuation does not play a dominant role [61]. The 
temperature makes a comparable contribution to microwave 
C? near the surface, particularly in the afternoon. For the 
acoustic and optical waves, the C? primarily depends on the 
temperature fluctuations. 


F Surface Roughness 


The wind profile in the surface layer can be described as 
[50], [51] 


Ux Z 
go I ~ Ym ; Lijo < 
r (n= w z/ 0 


a 


(29) 
Ux Z Z 
—— |i s m Le ’ Lro > 
Kay | - Zo y ( Dei ) z] j 


where z; 1s the height of boundary layer, and Ym 1s calculated 
as [50] 


a d(z/Lino 
a= (1 — bm) Am) ao) 


o/ Lane z/Lmo 


By substituting (23) into (30) and ignoring z9/Lmo, (30) can 
be expressed as 


—5z/Lmo, 2) ne 0 


In { [(1 + 29) /2] [(1+21)/27 
—2tan-12, + 7/2, 
Ym a 0 P Zi Dio > =2 (31) 
(3/2) In [(a5 + v2 + 1)/3] 
—/3tan~ |(2x%2 + 1)/V3] + 7/3, 


Z/ Imo < —2 


where xı = (1 — 16z/Lmo)!/4 and z2 = (1 — 15z/Lmo)!”?. 

A smooth sea surface can be approximated as a perfectly 
conducting surface, and a rough sea surface can be modeled 
with the Miller-Brown roughness reduction factor. The wind 
speed at 10 m height, U10, is used to represent the on-shore 
wind under a near-neutral atmospheric condition. 


G. Evaporation Duct Height 


Duct height variation with wind speed has been reported in 
[62]. The coastal areas are especially rich in super-refractive 
layers and ducts that affect microwave propagation [63]. In 
[64], cumulative distribution functions of duct height are 
derived in three locations. A larger variation in duct height is 
observed around the Lucinda data than in the Coral Sea. When 
sea breezes are the dominant wind flow, duct heights can be 
highly variable over a 24 hour cycle. The ducting events in 
the Coral Sea occurred under stable weather conditions, with 
little variation in wind speed; while the ducting events in the 
Western Pacific Ocean shows the greatest variation in duct 
height because of the highly variable weather. 
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Fig. 1. Modified refractivity profiles in an evaporation duct with M (0) = 
330, zq = 29 m, : very unstable, — — —: unstable, — —: near-neutral 
unstable, —.—: neutral, —o—: near-neutral stable, —«—: stable, —A—-: 


very stable, Lmo and zo of table I are used. 


IV. SIMULATION AND DISCUSSIONS 


Fig.l shows examples of modified refractivity profiles ex- 
pressed in (16), in an evaporation duct with zg = 29 m, under 
seven atmospheric conditions as listed in Table I. 

The difference of M values at the duct height and on the 
sea surface changes significantly under different atmospheric 
conditions. The atmospheric condition is well correlated to the 
Monin-Obukhov length, Ly, ., and the latter can be estimated 
using the Richardson number, Ff, as in (17). The Richardson 
number is a function of the Bowen ratio, as in (18), and the 
latter is expressed in (19). When the atmosphere becomes 
more unstable, the humidity gradient, Oq/Oz, increases in 
magnitude, leading to a stronger duct, or a larger difference 
of M values at the duct height and on the sea surface. 

A one-dimensional (vertical) Kolmogorov spectrum is 
adopted in [23], hence the phase fluctuation of waves prop- 
agating in the horizonal direction is not considered. This 
approximation is reasonable if the step size, Ax, in the SSF 
model is much larger than Lon. However, when the atmosphere 
becomes more unstable, the horizontal outer scale increases 
and may become comparable to Az [55]. In that case, a two- 
dimensional spectrum is more suitable. 

In this work, the two-dimensional anisotropic Kolmogorov 
spectrum of refractive-index fluctuation in (15) is adopted, 
with the outer scale of the turbulence in the surface layer 
expressed in (24), the structure constant in an evaporation duct 
expressed in (28), and the Loz/Loy ratio under z/Lmo < 0.5 
expressed in (26). 

Figs.2, 3 and 4 show the variations of the vertical outer 
scale, Loy, the structure constant, C?, and the Lop, /Loy ratio, 
respectively, covering all the seven atmospheric conditions 
listed in Table I. Since evaporation ducts are formed only 
tens of meters above the sea surface, the meteorological 
models which are valid in the surface layer are sufficient to 
characterize these ducts. 

The outer scales under unstable conditions are much larger 
than those under stable conditions, which has been confirmed 
by observations [55]. Under unstable conditions, the outer 
scale grows with height, as the air convection tends to reduce 
the wind shear. Under stable conditions, the outer scale tends 
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to decrease with height, as the stratification process tends to 
increase the wind shear. Accompanied by larger refractive- 
index fluctuation and larger outer scales, the structure constant 
under unstable conditions becomes much larger than that under 
stable conditions. Fig.1 implies that, under unstable conditions, 
stronger ducting effect allows the wave to propagate with a 
larger grazing angle, which in turns suffers stronger scattering 
effect by the rough sea surface [19]. Figs.2 and 3 indicate that, 
under unstable conditions, C2 and Lo become larger, hence the 
air turbulence gets stronger. The surface layer usually appears 
in the lowest 10 % of the boundary layer. At the height above 
zi/10, the height profiles of C2, Los and Lon/Lov, as shown 
in Figs.1, 2, 3 and 4, respectively, may not be accurate enough. 
Since the propagation properties in and near the surface layer 
are of major concern, less accurate profiles at higher altitudes 
should cause little problem. 


10° 








1 0 20 40 60 80 100 
Height (m) 
Fig. 2. Vertical outer scale Lo, over height, : very unstable, — — —: 
unstable, — —: near-neutral unstable, —.—: neutral, —-o—: near-neutral 


stable, —o—: stable, —A—: very stable, [mo of table I is used. 





0 20 40 60 80 100 
Height (m) 


Fig. 3. | Structure constant C2 over height, : very unstable, — — —: 
unstable, — -—: near-neutral unstable, —.—: neutral, —-o—: near-neutral 
stable, —o—: stable, —A—: very stable, Lmo and zo of table I are used. 





As shown in Fig.4, the Lo, /Loy ratio in a very stable air 
(z/Lmo => 0.5) becomes unpredictable when z > 0.5L no, yet 
its average value above 0.5Lmo 1s close to those in a slightly 
unstable air [50], [58]. In this work, the Lon/Loy ratio in 
z > 0.5L mo 18 approximated as that of a neutral air, namely, 
Lon/Lov = 1.592. The reference wind speed at 10 m above 
the sea is estimated with (29). 
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Fig. 4. Anisotropy ratio Lon /Lo» over height, : very unstable, — — —: 
unstable, — —: near-neutral unstable, —.—: neutral, —o—: near-neutral 
stable, —o—: stable, —A—: very stable, [mo of table I is used. 





Referring to [65], we choose zz = 29 m, % = 4 m, 
f = 10.6 GHz, in the subsequent simulations. The wave of 
horizontal polarization, which is more significantly affected 
than the vertical polarization in the two-dimensional case, 
is chosen. Radio signals in the 10.6 GHz band are often 
trapped in an evaporation duct above the sea surface, hence 
can be transmitted well beyond the horizon. The duct height 
of zg = 29 m 1s a typical number observed in the afternoon. 

Figs.5 to 8 show the path-loss over range and height in 
evaporation ducts under very unstable (vu), neutral (n), near- 
neutral stable (ns), and stable (s) atmospheric conditions, 
respectively, as listed in Table I. The frequency is f = 10.6 
GHz, the duct height is zg = 29 m, and the transmitter antenna 
height is z = 4 m. 

As shown in Fig.l, when the atmospheric condition pro- 
gresses from stable to very unstable, the slope, |dM/dz], 
within the duct increases significantly; and the slope outside 
the duct increases slightly. By comparing case (a) in Figs.5 
to 8, different characteristics in the path-loss distribution is 
observed between the unstable (Lmo < 0) and the stable 
(Lmo > 0) conditions. Although ducts appear under all 
atmospheric conditions, the contrast of path-loss within the 
duct and outside the duct is more obvious under unstable 
conditions. 

The path-loss distribution within the duct displays an inter- 
ference pattern between the propagating wave and the reflected 
wave from the smooth sea surface. The period of interference 
pattern within the duct under very unstable condition is shorter 
than those under the other three conditions, which is attributed 
to the larger grazing angle under the former condition. The 
interference pattern seems to undermine the ducting effect. 
Under unstable condictions, the effect of rough sea surface 
increases as the grazing angle, x, in (9), is increased. Mean- 
while, the turbulence effect is increased as (O(n) /Oz)? in (27) 
is increased. 

When the atmospheric condition is close to neutral (R; ~ 
0), (18) implies that the wind shear, |Ov;, /Oz|, becomes larger, 
hence the wind speed, Uj9, also becomes higher. Note that in 
(8) and (9), the effect of rough sea surface is dominated by 
the grazing angle, x, and the wind speed, Uio. 
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Fig. 5. Path-loss over range and height in a very unstable evaporation duct 
above rough sea surface, with z+ = 4 m, zą = 29 m, f = 10.6 GHz, 
Lmo = —73 m, us = 0.22 m/s, zo = 6.2 x 1075 m, and U10 = 6.4 m/s, 
(a) smooth sea surface, without turbulence, (b) rough sea surface, without 
turbulence, (c) rough sea surface, with turbulence. 


By comparing case (b) and case (a) in Figs.5 to 8, the path- 
loss near the sea surface at long ranges in case (b) (with rough 
sea surface) is always larger than that in the corresponding case 
(a) (with smooth sea surface). On the other hand, the path-loss 
at the same range but higher altitudes in case (b) becomes 
smaller than its counterpart in case (a). This phenomenon 
of energy leakage (from the duct to the air above) may be 
attributed to the scattering by the rough sea surface. The 
leakage phenomenon is more obvious under very unstable and 
neutral conditions. 
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Fig. 6. Path-loss over range and height in a neutral evaporation duct above 
rough sea surface, with z = 4m, z = 29m, f = 10.6 GHz, Lmo = 
—1531 m, ux = 0.4 m/s, z = 19.6 x 1075 m, and Uio = 10.81 m/s, 
(a) smooth sea surface, without turbulence, (b) rough sea surface, without 
turbulence, (c) rough sea surface, with turbulence. 


Eqns.(8) and (9) suggest that the roughness effect of sea 
surface becomes more obvious as the grazing angle increases 
and/or the wind speed increases. As shown in Fig.1, a larger 
gradient of refractive index is observed near the sea surface 
when the atmosphere is more unstable, causing a larger grazing 
angle. Based on (29), the wind speed, Uj9, becomes larger 
when the atmospheric condition is close to neutral. 

The path-loss distributions under unstable (u) and near- 
neutral unstable (nu) conditions appear similar to those un- 
der very unstable (vu) condition and neutral (n) condition, 
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Fig. 7. Path-loss over range and height in a near-neutral stable evaporation 
duct above rough sea surface, with z = 4 m, zq = 29 m, f = 10.6 GHz, 
Lmo = 314 m, ux = 0.23 m/s, zo = 6.3 x 1075 m, and U10 = 6.98 m/s, 
(a) smooth sea surface, without turbulence, (b) rough sea surface, without 
turbulence, (c) rough sea surface, with turbulence. 


respectively, hence are not presented. Under very stable (vs) 
condition, neither the rough sea surface nor the air turbulence 
generates noticeable difference in the path-loss. Hence, it is 
not presented, either. 


As shown in Fig.3, the structure constant decreases with 
height. This tendency is predictable from (27), using the 
profiles of |dM/dz| and Loy shown in Figs.1 and 2, re- 
spectively. Larger value of C? implies stronger turbulence 
intensity, especially under neutral and unstable conditions. 
Although larger value of C? is observed near the sea sur- 
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Fig. 8. Path-loss over range and height in a stable evaporation duct above 
rough sea surface, with z = 4 m, zy = 29 m, f = 10.6 GHz, Lmo = 85 
m, ux = 0.15 m/s, zo = 2.9 x 107° m, and Ujg = 5 m/s, (a) smooth 
sea surface, without turbulence, (b) rough sea surface, without turbulence, (c) 
rough sea surface, with turbulence. 


face, under all conditions; the turbulence effect appears more 
obvious at higher altitudes attributed to the outer scales in the 
Kolmogorov spectrum of (15). As the atmospheric condition 
progresses from neutral to very unstable, the turbulence creates 
strong interference in the upper right quarter of the path-loss 
distribution, and significantly increases the path-loss in that 
area. The path-loss in the duct is increased by the order of 10 
dB at the range of 100 km. 


Fig.9 shows the standard deviation of path-loss below the 
duct height of zg = 29 m, based on the data presented in Figs.5 
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Fig. 9. Standard deviation of path-loss within the duct (zg = 29 m), 
effect of rough sea surface, — — —: effect of turbulence, (a) very unstable 
condition (vu), (b) neutral condition (n), (c) near-neutral stable condition (ns), 
(d) stable condition (s). 





to 8. The standard deviation of path-loss, at a given range, 
attributed to the rough sea surface is derived by taking the root- 
mean-square (rms) difference between the data with smooth 
sea surface (case (a)) and those with rough sea surface (case 
(b)). Similarly, the standard deviation of path-loss, at a given 
range, attributed to the turbulence is derived by taking the 


10 


root-mean-square difference between the data with rough sea 
surface without turbulence (case (b)) and those with turbulence 
(case (C)). 

Strong rough sea surface effect is observed under very 
unstable and neutral conditions, attributed to the strong ducting 
effect under very unstable condition, as shown in Fig.l, and 
the high wind speed under neutral condition, in (17) and 
(18). The rough sea surface contributes no obvious effect 
under near-neutral stable and stable conditions, due to the low 
wind speed and weak ducting. As the atmospheric condition 
progresses from stable to very unstable, additional path-loss 
due to turbulence gradually increases, as can be predicted from 
Figs.2 and 3. 

The main contribution of this work is to simulate the 
combined effects of rough surface and inhomogeneous atmo- 
sphere to the propagation in evaporation ducts. Although many 
experimental works have been successfully conducted before, 
it still requires some theoretical analysis and simulation to 
differentiate the separate effects of the aforementioned mech- 
anisms. The parameters chosen in this work are closely tied 
to the previous experiments, hence the simulation results are 
complementary to those works. A two-dimensional simulation 
scenario takes less computational burden to study these effects. 
The simulation approach in this work can be extended to three- 
dimensional cases (with two-dimensional rough surfaces), only 
at a heavier computational load. 

The paraxial approximation leads to more accurate pre- 
dictions at high frequencies and in the forward direction. 
In this work, the frequency is chosen to be 10.6 GHz, and 
the propagation direction in evaporation ducts is pretty much 
forward. When using the split-step Fourier method, the step 
size in the propagation direction is restricted by the minimum 
scale length of the inhomogeneous atmosphere. The step size 
becomes very small if the tiny details of the atmosphere is 
to be included, then the computational time will be increased 
and numerical instability may occur at long distances. 


V. CONCLUSION 


A proper set of physical and mathematical models have 
been compiled to analyze and simulate wave propagation in 
evaporation ducts, under the scattering effect of rough sea 
surface and the air turbulence. All the relevant parameters are 
clearly connected to seven possible atmospheric conditions. 
A systematic approach, including the range of parameters, 
has been proposed to simulate wave propagation in close- 
to-realistic evaporation ducts. Simulations have also been 
conducted to demonstrate the effectiveness of this model and 
approach. 

The split-step Fourier propagation algorithm has been ap- 
plied to analyze the propagation properties in evaporation 
ducts above a rough sea surface, incorporating the turbulence 
effect. The PJ model of modified refractivity under different 
atmospheric conditions is adopted in our model, which renders 
significant difference of path-loss distribution under different 
atmospheric conditions. 

The scattering effect of the rough sea surface is simulated 
by introducing a correction factor of surface roughness into the 
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boundary condition. When the atmosphere is more unstable, 
the interference pattern within a duct bears a shorter period 
and a wider range of path-loss variation. 

The turbulence is described with a two-dimensional Kol- 
mogorov spectrum, in which the structure constant of refrac- 
tive index fluctuation, the outer scales of turbulence, and the 
anisotropy of the outer scales are well connected to different 
atmospheric conditions. 

Significant effects of rough sea surface and air turbulence 
are observed under neutral (|Lmo| ~ 00) and unstable (Lino < 
0) conditions. The effect of rough sea surface is less significant 
under stable (Lm. > 0) conditions, accompanied by low 
grazing angle and low wind speed. The effect of air turbulence 
gradually increases as the atmospheric condition progresses 
from stable to unstable, accompanied by larger eddy size and 
stronger turbulence intensity. 
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Editorial Comment 


In this work the authors have addressed a very challenging problem, namely that of evaluating 
the effect of turbulence on wave propagation in evaporation ducts above a rough surface. Several 
different available models have been used to simulate the problem, and the results of these 


simulations are presented as well as discussed. 


